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Appendix A 



<?xml version= ,, 1.0 n encoding = "UTF-8" ?> 

- <model name= ,, FieldKorosNoyesModel"> 

- <notes> 

5 <hl>Field-Koros-Noyes Model of BZ Reaction</hl> 

- <table border= ,, 0" cellspacing= ,, 0 M cellpadding = "2 ,, > 

- <thead> 

- <tr> 

<th align="left" valign= "middle" bgcolor="#eeeeee M >Citation</th> 

10 </tr> 

</thead> 

- <tbody> 

- <tr> 

- <td> 

15 R.J. Field and R.M.Noyes,J.Chem.Phys.60, 1877 (1974) ; 

R.J. Field, E.Koros,R-M.Noyes,JACS 94,8649 (1972);R.J. Field, R.M.Noyes,Nature 
237,390 (1972) This implementation is taken manufactured by J.D. Murray, 
"Mathematical Biology" (1989) page 181. 

<a href= ,,M /> 
20 </td> 
</tr> 
</tbody> 
</table> 

- <table border="0 M cellspacing="0" cellpadding="2 ,, > 
25 <thead> 

- <tr> 

<th align="left" valign = "middle" bgcolor="#eeeeee M >Description</th> 

</tr> 
</thead> 
30 <tbody> 

- <tr> 

<td> Field Noyes Version of Belousov- Zhabotinsky Reaction. Br03 is held 

constant; HOBr is typically ignored, and can be replaced by an empty- set. The 
stoichiometry f is typically taken as 1/ 2 or 1 (denominator 1 or 2 in SBML) 

35 .</td> 
</tr> 
</tbody> 
</table> 

- <table border= ,, 0" cellspacing = M 0" cellpadding = ,, 2 M > 
40 <thead> 

- <tr> 

<th align = N left" vaiign="middle" bgcoIor="#eeeeee">Rate constant </th> 
<th align = M left" valign = n middle" bgcolor="#eeeeee">Reaction</th> 

</tr> 
45 </thead> 

- <tbody> 

- <tr> 

<td>kl = 1.3</td> 

<td>Br + Br03 -> HBr02 + HOBr</td> 

50 </tr> 
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- <tr> 

<td>k2 = 2000000</td> 

<td>Br + HBr02 -> HOBr*2</td> 

</tr> 
§ <tr> 
<td>k3 = 34</td> 

<td>Br03 + HBrOZ -> Ce^2 + HBr02 yv 2</td> 

</tr> 

- <tr> 

10 <td>k4 = 3000</td> 

<td>HBr02*2 -> Br03 + HOBr</td> 

</tr> 

- <tr> 

<td>k5 = 0.02</td> 
15 <td>Ce -> BrM</td> 
</tr> 
</tbody> 
</table> 

- <table border= ,, 0 H cellspacing = "O n cellpadding="2 n > 

29 <thead> 

- <tr> 

<th align="left" valign = "middle" bgcolor="#eeeeee">Variable</th> 
<th align='ieft" valign = "middle" bgcolor="#eeeeee">IC </th> 
<th align="left" valign = ,, middle" bgcolor= n #eeeeee">ODE</th> 

25 </tr> 

</thead> 

- <tbody> 

- <tr> 
<td>Br</td> 

30 <td>0.003</td> 

<td>Br'[t] == -(kl*Br[t]*Br03[t]) + f*k5*Ce[t] - k2*Br[t]*HBr02[t]</td> 

</tr> 

- <tr> 
<td>Ce</td> 

35 <td>0.05</td> 

<td>Ce'[t] == -(k5*Ce[t]) + 2*k3*Br03[t]*HBr02[t]</td> 
</tr> 

- <tr> 
<td>HBr02</td> 

40 <td>0.00K/td> 

<td>HBr02'[t] == kl*Br[t]*Br03[t] - k2*Br[t]*HBr02[t] + 
k3*Br03[t]*HBr02[t] - k4*HBr02[t] /v 2</td> 
</tr> 

- <tr> 

45 <td>HOBr</td> 
<td>0</td> 

<td>HOBr'[t] == kl*Br[t]*Br03[t] + 2*k2*Br[t]*HBr02[t] 
+ k4*HBr02[t]^2</td> 

</tr> 
50 </tbody> 
</table> 
</body> 
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</notes> 

- <MstOfCompartments> 
<compartment name="BZ" /> 

</MstOfCompartments> 
5 <listOfSpecies> 

<specie name="Br" initialAmount="0.003" compartment="BZ" boundaryCondition="false" 

/> 

<specie name="Br03" initialAmount= ,, 0-l" compartment="BZ" boundaryCondition = "true" 
/> 

10 <specie name="Ce" initialAmount="0.05 n compartment="BZ" boundaryCondition = "false" 
/> 

<specie name= ,, HBr02 ,, initialAmount= M 0.001" compartment="BZ M 

boundaryCondition="false" /> 
<specie name="HOBr" initialAmount= n O" compartment= ,, BZ M boundaryCondition = "false" 
15 /> 

</listOfSpecies> 

- <NstOfReactions> 

- <reaction name="Reactionl" reversible="false"> 

- <MstOfReactants> 

20 <specieReference specie="Br" /> 
<specieReference specie="Br03" /> 
< / 1 i stOf Rea eta n ts > 

- <MstOfProducts> 
<specieReference specie="HBr02 M /> 

25 <specieReference specie="HOBr" /> 
</IistOfProducts> 

- <kineticLaw formufa="Br*Br03*kl ,, > 

- <listOfParameters> 

<parameter name="kl" value="l-3" /> 
30 </listOfParameters> 
</kineticLaw> 
</reaction> 

- <reaction name="Reaction2" reversible= ,, false"> 

- <listOfReactants> 

35 <specieReference specie="Br" /> 
<specieReference specie="HBr02" /> 
</listOfReactants> 

- <listOfProducts> 

<specieReference specie= ,, HOBr" stoichiometry= M 2 M /> 
40 < / 1 i s tOf P ro d u ct s > 

- <kineticLaw formula= ,, Br*HBr02*k2"> 

- <listOfParameters> 

<parameter name="k2 M value="2000000" /> 
</l istOf Pa ra meters> 
45 </kineticLaw> 
</reaction> 

- <reaction name= , 'Reaction3 ,, reversible="false"> 

- <MstOfReactants> 
<specieReference specie= M Br03" /> 

50 <specieReference specie="HBr02 M /> 
</listOfReactants> 

- <MstOfProducts> 



3 



Docket No.: MWS-108 

<specieReference specie="Ce" stoichiometry="2" /> 
<specieReference specie= M IHlBr02 ,, stoichiometry= ,, 2" /> 
</MstOfProducts> 

- <kineticLaw formula = ,, Br03*HBra2*k3"> 
5 <MstOfParameters> 

<parameter name="k3" value= ,, 34 n /> 
</listOfParameters> 
</kineticLaw> 
</reaction> 

19 <reaction name = M lReacttioni4" reversible= ,, 1False ,, > 

- <MstOfReactants> 
<specieReference specie= ,, IHIBr02 ,, /> 

</MstOfReactants> 

- <listOfProducts> 

15 <specieReference specie="Br03" /> 
<specieReference specie^'HOBr" /> 
</listOfProducts> 

- <kineticLaw formula = ' , HBr02*2*k4 H > 

- <MstOfParameters> 

20 <parameter name="k4" value="3O(Q>0" /> 

</l istOf Pa ra m ete rs > 

</kineticLaw> 

</reaction> 

- <reaction name = ,, Reaction5 n reversible="ffalse"> 
25 <listOfReactants> 

<specieReference specie="Ce" /> 
</listOfReactants> 

- <listOfProducts> 

<specieReference specie="Br" stoichiometry="l M denominator="2" /> 
30 </listOfProducts> 

- <kineticLaw formula= M Ce*k5"> 

- <listOfParameters> 

<parameter name="k5" value="0.02" /> 
<parameter name="1F" value^'O.B" /> 
35 </MstOfParameters> 
</kineticLaw> 
</reaction> 
</MstOfReactions> 
</model> 
40 </sbml> 
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